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Large networks of integrate-and-fire (IF) model neurons 
are often used to simulate and study the behaviour of bio- 
logically realistic networks. However, to fully study the 
large network behaviour requires an exploration of large 
regions of a multidimensional parameter space. Such 
exploration is generally not feasible with large network 
models, due to the computational time required to simu- 
late a network with biologically significant size. To circum- 
vent these difficulties we use a mean-field approach, based 
on the work of [1], 

We consider a sparsely coupled, excitatory network of 
10,000 Izhikevich model neurons [2], with Destexhe- 
type synapses [3]. The cellular models were fit to hippo- 



campal CA1 pyramidal neurons and have heterogeneous 
applied currents with a normal distribution. We derived 
a mean-field system for the network which consists of 
differential equations for the mean of the adaptation 
current and the synaptic conductance. 

As CA1 is an area that displays prominent theta oscilla- 
tions [4], we used the mean-field system to study how the 
frequency of bursting depends on various model para- 
meters. Figure 1A shows an example study. These studies 
were successful in guiding numerical simulations of the 
large network. When parameter values determined from 
the mean-field analysis are used in a large network simula- 
tion, bursting of the predicted frequency occurs (Figure IB). 
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Figure 1 A. Mean-field prediction of the bursting frequency as a function of the unitary synaptic conductance and the mean applied current. 
B. Simulation of a network of 10,000 neurons with unitary conductance 0.058 nS and mean current 250 pA, showing an oscillation in the theta 
frequency range as predicted by the mean-field system of equations. 
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